Water temperature dependence of single bubble sonoluminescence 



Sascha Hilgenfeldt 1 , Detlef Lohse 1 , and William C. Moss 2 
1 Fachbereich Physik der Universitdt Marburg, Renthof 6, D-35032 Marburg 
2 Lawrence Livermore National Laboratory, Livermore, CA 94550 
(Phys. Rev. Lett., February 16, 1998) 



The strong dependence of the intensity of single bubble 
sonoluminescence (SBSL) on water temperature observed in 
experiment can be accounted for by the temperature depen- 
dence of the material constants of water, most essentially of 
the viscosity, of the argon solubility in water, and of the va- 
por pressure. The strong increase of light emission at low wa- 
ter temperatures is due to the possibility of applying higher 
driving pressures, caused by increased bubble stability. The 
presented calculations combine the Rayleigh-Plesset equation 
based hydrodynamical/chemical approach to SBSL and full 
gas dynamical calculations of the bubble's interior. 



One of the remarkable features of single bubble sono- 
luminescence (SBSL) jy,g] is the sensitivity of the light 
emission to the water temperature experimentally found 
by the UCLA group c£ figure [l]. To obtain these 

results, Barber et al. proceeded as follows (refs. J2|, [| and 
B. Barber, private communication, August 1997): Wa- 
ter was cooled to a temperature of 2.5°C and completely 
degassed. An air pressure overhead of 150 Torr, corre- 
sponding to about 20% of gas saturation, was adjusted 
and SL experiments were performed, still at 2.5°C. Then 
the water was heated to 20°C without readjusting the gas 
concentration, and the SL experiment was repeated. Fi- 
nally, the same measurement was performed after heat- 
ing the water to 33°C. At all three temperatures, the 
forcing pressure amplitude P a of the driving sound field 
was adjusted in order to give maximum light intensity, 
while maintaining bubble stability against fragmentation 
(stable SL). According to the "waterfall plots" shown e.g. 
in ref. Q , the highest light intensity always corresponds 
to the largest achievable driving pressure (P™ ax ). The 
experimental values for p™ ax are shown in figure g, 

The bubble radius dynamics R(t) was detected in these 
experiments (using Mie scattering techniques) and fit- 
ted by the Rayleigh-Plesset (RP) equation [||. In this 
fit, several parameters were allowed to vary: the driving 
pressure amplitude P a , the ambient radius i?o (radius of 
the bubble under normal conditions) and also the surface 
tension a of the air/ water interface and the water viscos- 
ity v\ |2|, H . We consider the treatment of the latter two 
quantities as free parameters arbitrary, as a and v\ are 
well defined material constants, but are varied in [pL||] 
by as much as a factor of four and beyond. Fitting a 
and vi may be avoided if a more realistic model for the 
thermal behavior of the gas inside the bubble is applied: 
in the RP equation, the internal gas pressure is taken to 
vary polytropically with volume, with an effective expo- 
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FIG. 1. Maximally achievable number N max of SBSL 
photons for given water temperature. The experimental data 
(open boxes) are taken from figure 1 of The theoretical 
data (filled circles) result from the full gas dynamical calcu- 
lations performed with the algorithm described in 

nent 7. In [||^] 7 changes abruptly from its isothermal 
to its adiabatic value at R(t) — Rq. Taking instead a 
more realistic model 7=1 throughout the whole oscilla- 
tion with the exception of the immediate vicinity of the 
strong bubble collapse results in a satisfactory fit with 
the physical values for a and v\ . 

Because of the complications in the fits of refs. j|, [|, 
the resulting data should be read with some care. This is 
also reflected in figure |^, where we display the data from 
these two references for the expansion ratio (maximum 
radius R m ax divided by Rq): they show large deviations 
at otherwise unchanged parameters. The expansion ra- 
tio is a quantity closely related to the violence of collapse 
and therefore, presumably, to the intensity of energy con- 
centration and light emission Q . It is therefore puzzling 
that the same light intensity has been observed in refs. 
[^],|3| in spite of the different expansion ratios reported in 
figure [| 

The central claim of this paper is that the observed 
dependence on water temperature T in figures [l] - [| can 
be accounted for by the T dependence of the material 
constants of water which are listed in table |[ In our 
analysis the most relevant effects are the temperature 
dependence of the gas solubilities (as already conjectured 
in ref. Q) and of the viscosity; also, the temperature 
dependence of water vapor contributes. The variations 
of these material constants with temperature conspire to 
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FIG. 2. Maximum forcing pressure for stable SBSL as a 
function of the water temperature. The experimental data 
(open boxes) are from figure 1 of ]j| . The filled circles are 
the theoretical values following from figure ^j. 
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TABLE I. Material constants of water as a function of 
temperature ||. From top to bottom: surface tension a, den- 
sity p, kinematic viscosity vi, vapor pressure p vap , speed of 
sound a, and solubilities of argon and nitrogen in water Cq t , 
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allow for larger shape stable bubbles and larger forcing 
at lower temperatures, resulting in more light emission. 

When dissolving air with its various constituents in 
water, only the argon concentration c^ r (far away from 
the bubble) is relevant in the SL regime, as SL air bub- 
bles rectify argon 0. The degree of saturation c^/cq 
is one of the central parameters to determine the diffu- 
sive equilibrium (ambient) radius of SL bubbles. This 
equilibrium radius can be calculated within the hydro- 
dynamical/chemical approach to SL j7j-S], which takes 
into account the mass exchange between the bubble in- 
terior and the liquid due to (i) diffusion and (ii) disso- 
ciation of molecular gases, which leads to argon rectifi- 
cation in the bubble. Here we follow that approach and 
calculate the phase diagrams in Rq — P a space, result- 
ing from the diffusional/chemical stability of the bubble 
for varying temperature, i.e., varying material constants 
and also varying c^ r /c^ r (T). The calculations (based on 
Rayleigh-Plesset dynamics) are described in detail in refs. 



FIG. 3. Experimental data for the expansion ratio 
Rmax/Ro at maximal ambient radius R<? ax as a function of 
the water temperature taken from figure 1 of || (open trian- 
gles) and figure 47 of || (open boxes). Also shown are the 
theoretical values following within the presented theoretical 
approach (filled circles). The lower part of the figure shows 
the maximal pressure achieved in the bubble within the hy- 
drodynamical/chemical approach. 



|7j . We choose the material constants given in table [| the 
frequency 26.5kHz as in experiment [0, and an ambient 
pressure of 1 atm. We modeled air as a mixture of 99% 
nitrogen and 1% argon (in the gas phase above the liq- 
uid). At a given total gas concentration of Coq/cq = 20% 
for T — 2.5°C, the corresponding argon concentrations 
for the three temperatures are readily calculated from ta- 
ble |: c^/4 r (2.5°C) = 0.20%, /c Ar {20°C) = 0.29%, 
and c^7c^ r (33°C) = 0.36%, respectively. 

The resulting equilibrium curves in phase space are 
shown in figure [|. Diffusively stable bubbles are possible 
on the branches B and C 0; however, according to the 
(qualitative) energy focusing condition \R\/c g > 1 
HO] (where c g is the speed of sound in the gas bubble) 
most bubbles on branch B will not be able to emit SL 
light, in agreement with recent experiments 
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A key issue is that the upper limit in both P a and i?o 
on branch C is given by the parametric shape instabil 
ity of the bubble wall 



J)|,y,|l4|. The shape instability 
depends on the viscosity of the water [jl3| , which strongly 
increases with decreasing water temperature, thus stabi- 
lizing the bubble. For the values given in table || we find 
(within the approximations of ||]) that shape instabil- 
ity sets in at an ambient radius of roughly Rq = 3.5/im, 
Afj,m, and 5/j,m for 33°C, 20°C, and 2.5°C, respectively, 
see the dashed lines in figure |[ 

From figure |] one immediately realizes that the bubble 
can be driven harder (larger p™ aa: ) for lower tempera- 
tures. The values P™ ax from figure [| are displayed in 
figure ||, showing good agreement with experiment. 

On decreasing the water temperature, the material con- 
stants of water are changed according to table R this has 
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FIG. 4. Phase diagram in the Ro — P a parameter space. 
The solid lines show the stable diffusive equilibria for three 
different temperatures. The bubbles along the equilibria C 
contain essentially pure argon. The bubbles along the equi- 
libria B have a slightly enhanced argon concentration - these 
bubbles are in equilibrium because the mass gain by rectified 
diffusion and the mass loss by chemical reactions balance Q. 
Left of each curve C bubbles are shrinking, right of each curve 
C they are growing, as indicated by the arrows. - The dashed 
lines (shape inst.) show where the bubbles become spherically 
unstable because of the parametric shape instability Q . The 
intersection with the diffusive stability line defines p™ aa: and 
R^ ax (circles). 



the following effect on the curves in figure § (i) The 
water viscosity v\ is increased; as stated above, vi deter- 
mines the damping of bubble shape oscillations; there- 
fore, bubbles are more stable; the stability line shifts to- 
wards larger Rq. (ii) The relative argon concentration 
c oo I c q t (T) becomes smaller; in order to counteract dis- 
solution, the bubbles must oscillate more violently to rec- 
tify gases. Thus, diffusive equilibria are reached only at 
higher P a , and curves B and C are shifted to the right, 
allowing for larger p™ ax . (hi) The water vapor pressure 
Pvap decreases; this, too, slightly shifts curves B and C to 
the right. The reason is that the reduced total pressure 
inside the bubble leads to smaller R max - The resulting 
reduced gas rectification has to be compensated by larger 
P a . (iv) The small temperature dependences of a, p, and 
ci have hardly any effect. 

In figure |3J we include the theoretical expansion ra- 
tio Rmax/ Ro\ii™ a * calculated with the Rayleigh-Plesset 
equation for p™ ax and R™ ax determined from figure 4. 
Here, the agreement with the data from the fits of 0, 2] 
is not as good. We speculate that part of the discrepancy 
originates from the fit procedure (to obtain Rq) of refs. 
[^],|3| (sec our criticism above) and suggest direct mea- 
surements of R m ax/ Ro, either by Mie scattering at two 
different frequencies or by employing the method devel- 



oped by Holt and Gaitan |ll| . 

One may wonder why a relatively modest change in 
P a and R Q nax (cf. figs. [| and ||) leads to considerably en- 
hanced light emission, cf. fig. @. The reason is twofold: 
(i) With increasing P a the bubble collapse becomes much 
more violent and therefore the gas inside the bubble un- 
dergoes stronger heating. The violence of the bubble 
collapse can be quantified by maximal pressure p max 
reached inside the bubble (within the hydrodynamical 
approach), cf. fig. ^, lower part. The observed increase 
with decreasing water temperature is a direct conse- 
quence of the observed increase in R ma x/Ro (upper part 
of fig. |3|): For bubbles with a larger expansion ratio, 
more potential energy can be converted into kinetic en- 
ergy during the collapse, compressing the bubble more 
strongly, (ii) Larger R™ ax means that a larger number 
of gas molecules oc (R T nax ) 3 can be heated, which also 
leads to more light emission. 

To complete our analysis, we would like to calculate 
the light intensity and compare it with the experimental 
data in figure ^. The light intensity can of course not be 
obtained from the RP approach to SBSL. It sensitively 
depends on the exact temperature achieved in the bubble 
at the collapse which can only be estimated within the 
Rayleigh-Plesset bubble approach, using one of a number 
of simple models, e.g., adiabatic heating at the collapse. 

However, numerical codes that calculate the gas dy- 
namics of the bubble's interior and the resulting optical 
emission have been developed |4[|l5|. We use the model 
recently proposed by Moss et al. || to calculate the gas 
dynamics inside the bubble at the collapse and the result- 
ing optical emission. The model solves the gas dynamic 
equations for the conservation of mass, momentum, and 
energy in the gas bubble and the surrounding liquid, us- 
ing a finite differencing scheme for the assumed spherical 
growth and collapse of the bubble. Accurate high pres- 
sure/temperature equations of state are used to describe 
the water and the argon. Energy loss by thermal con- 
duction in the partially ionized gas created during the 
collapse is included in the model. Although light emis- 
sion is the primary diagnostic of SBSL, the energy loss 
by the radiation of light is negligible compared to that 
by thermal conduction. This simplifies the calculation of 
the emitted optical power and spectra, which are calcu- 
lated using the opacity and the calculated densities and 
electron temperatures of the partially ionized gas. Full 
details are given in [Q. The model was able to account 
for many experimentally found features of SBSL, includ- 
ing spectra and the sensitivity of the spectral intensity 
to the applied acoustic pressure (or, equivalently, R max )- 
This last result establishes the relevance of the model to 
the current analysis. We use the parameter pairs (p™ aa; ; 
R Q nax ) or equivalently {R m ax\R™ asc , i?™ ax ) resulting from 
above hydrodynamical/chemical approach for the three 
different analyzed water temperatures as input parame- 
ters for the gas dynamical code. The number of photons 
obtained in the wavelength window [20071771/75077,771] is 
compared to the experimental data in figure |l|, showing 
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reasonable agreement. 

To summarize, the experimentally found water tem- 
perature dependence of SBSL can consistently be 
accounted for by the water temperature dependence of 
the material constants of water: At lower water temper- 
ature bubbles are more stable and can therefore be forced 
more strongly, resulting in more light. Combined with a 
gas dynamical simulation of the bubble's interior, this 
approach can reproduce the three orders of magnitude 
increase in the observed photon number when changing 
the temperature from ~ 30°C to ~ 0°C. 

We suggest several control experiments to study the 
temperature dependence of SBSL. In all cases the exper- 
iments should be done with stable SBSL, whose regime in 
phase space can always (at least approximately) be cal- 
culated a priori, following the hydrodynamical/chemical 
approach. 

(i) The bubble dynamics R(t) should be measured as 
a function of the water temperature - this should re- 
veal that the water temperature dependence of SBSL is 
mainly a bubble dynamical effect. In particular, as al- 
ready mentioned above, the experimental data of figure 
|3| should be remeasured. 

(ii) In addition, as suggested in a recent preprint by 
Vuong, Fyrillas, and Szeri jl6), one should study how 
the bubble dynamics and the light intensity change (for 
fixed forcing pressure) with water temperature for xenon 
doped and helium doped nitrogen bubbles. The two no- 
ble gases should show different behavior, because c^ e 
strongly depends on the water temperature T, whereas 
c^ e only shows a weak dependence, see figure 3 of ref. 
Jl6| . This kind of experiment will be a useful diagnostics 
of the identity of the gas in the bubble [Q. It will also 
help to distinguish between the gas solubility effects on 
the light intensity and other bubble dynamical effects, 
e.g., those caused by the temperature dependence of the 
water viscosity and of the vapor pressure. 

(iii) Finally, instead of preparing the gas concentra- 
tion at one temperature (resulting in a water temperature 
dependent relative concentration c^ r /c^ r (T)), one could 
keep c^/cQ r constant at various temperatures, e.g. by 
using noble gas doped oxygen bubbles and controlling 
the actual gas content with oxyometry, cf. Gompf et al. 
fl7f . If in addition the forcing pressure amplitude is kept 
constant, the location of the shape instability and its de- 
pendence on water viscosity will also become irrelevant. 
Altogether, we therefore predict that the water tempera- 
ture dependence of SBSL light emission at otherwise fixed 
parameters P a , c^ff /cq t will be only weak. 
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